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ABSTRACT 

We present the characterization and initial results from the QUEST-La Silla AGN variability survey. 

This is an effort to obtain well sampled optical light curves in extragalactic fields with unique multi¬ 
wavelength observations. We present photometry obtained from 2010 to 2012 in the XMM-COSMOS 
field, which was observed over 150 nights using the QUEST camera on the ESO-Schmidt telescope. 

The survey uses a broadband filter, the Q-band, similar to the union of the g and the r filters, achieving 
an intrinsic photometric dispersion of 0.05 mag, and a systematic error of 0.05 mag in the zero-point. 

Since some detectors of the camera show significant non-linearity, we use a linear correlation to fit the 
zero-points as a function of the instrumental magnitudes, thus obtaining a good correction to the non¬ 
linear behavior of these detectors. We obtain good photometry to an equivalent limiting magnitude 
of r ~ 20.5. The astrometry has an internal precision of ~ 0'/l and an overall accuracy of 0"2 when 
compared to SDSS. Studying the optical variability of X-ray detected sources in the XMM-COSMOS 
field, we find that the survey is ~ 75 — 80% complete to magnitudes r ~ 20, and ~ 67% complete to a 
magnitude r ~ 21. Additionally, broad-line (BL) AGN have larger variation amplitude that non-broad 
line (NL) AGN, with ~ 80% of the BL AGN classified as variable, whilst only ~ 21% of the NL AGN 
resulted classified as variable. We also find that ~ 22% of objects classified as galaxies (GALs), are 
also variable. The determination and parameterization of the structure function ( SF norrn (r ) = At 1 ) 
of the variable sources shows that most BL AGN are characterized by A > 0.1 and 7 > 0.025. It is 
further shown that NL AGN and GAL sources occupying the same parameter space in A and 7 are 
very likely to correspond to obscured or low luminosity AGN. Our samples are, however, small, and 
we expect to revisit these results using larger samples with longer light curves obtained as part of our 
ongoing survey. 

Subject headings: Astrometry — Photometric techniques — quasars: general 


1. INTRODUCTION 

Despite the fact that variability is one of the defining 
characteristics of Active Galactic Nuclei (AGNs) we do 
not completely understand the mechanisms that drive 
such variation. Our understanding misses significant de¬ 
tails of how AGN variability at different wavelengths 
is rel ated (s ee e.g., lArevalo et al.l 120081 l2009t iLira et ahl 
[20Lfl. [20151) . and how physical properties of the cen¬ 
tral engine (e.g., luminosity, black hole mass, hardness 
ratio, optical colors, etc) are related to well defined 
variability properties of the system (e.g., characteristic 
timescale, variability amplitude, etc). From a cosmo¬ 
logical perspective, the strong evolution of AGNs opens 
the possibility to observe changes of the structure and 
feedin g of AGNs over cosmic time (see iShemmer et al.1 
na, and references therein). In fact, controversy re¬ 
mains on the redshift dependence of AGN variability 
since the observed wavelength and the minimum lu- 
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minosity i n a sample are highly c o rrelated with red- 
shift (see lYanden Berk et al.l 120041: Kelly et al.l 120091 : 
iMacLeod et al.ll2010HMorganson et al.ll2014 . and refer¬ 
ences therein). Finally, the property of being highly vari¬ 
able makes AGN selection by variability a very promising 
tool to find them. 

In order to study in detail the variability proper¬ 
ties of individual AGNs and to use this information to 
rule out different variability models to probe a wide 
range of timescales is required. Hence, long and in¬ 
tense campaigns are crucial. Fortunately, in recent years, 
surveys covering a significant part of the sky, revisit¬ 
ing the same regions on timescales from days to years, 
and containing a large sample of serendipitous objects 
blind surveys- are now becoming available as prede¬ 
cessors o f th e Lar ge Synoptic Survey Telescope (LSST; 
llvezic et al.l I2008T) . LSST will revisit each part of the 
of the southern sky (<5j2ooo <+34.5) approximately ev¬ 
ery three nights over 10 years, and will observe in six 
bands ( ugrizy ) to a limiting magnitude of r ~ 24.5 
mag. Among the LSST pr edecessors we ha ve: Sloan 
Digital Sky Survey (SDSS- dYork et ahl l2000h ) equato¬ 
rial Stripe-82 which covers roughly 280 deg 2 in the 
2Q h : 34 m -04 h : 08 m ri ght ascension range, and dec¬ 
linati on |< 5 j 2 ooo|< 1 - 266 ° (flvezic et al.l 12004 iSesar et al.l 
2007). This region was observed in five bands ( ugriz ) 
once to three times a year from 2000 to 2005 (SDSS- 
I), and then from 2005 to 2008 at an increased ca- 
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dence of 10 to 20 ti mes per year as par t of the SDSS- 
II supernova survey dFrieman et alJl2008f). The Catalina 
Real-Time Transien t Survey ('CRTS' TDrake et all 120091 
iGraham et al.ll2014lf which covers between —75° << 5 j 2 ooo 
<75°, observing in the V filter to a limiting magnitude 
of ~ 19 to 21.5 (depending on the telescope used). CRTS 
covers a total of 33, 000 deg 2 observing up to ~ 2500 deg 2 
per night, with 4 exposures per visit, separated by 10 
min, over 21 nights per lunation (i.e., re-visiting a field 
every 10 to 15 days). The Palomar T ransient Factory 
(PTF; IR.au et all[20091: ILaw et al.ll2009lf covers between 
— 15° < 6 j 2 ooo <85° using two broad-bandpasses, namely 
the Mould R band and g band, to a limiting magnitude 
of R ~ 21.0 mag and g ~ 21.6 mag with a cadence 
of 3 to 5 days. The Panoramic Survey T elescope and 
Rapid Response Syste m (PanSTARRS; iHodapp et al.1 
1200-11 [Tbiirv et ahl[ 20 12 ) covers all the northern sky down 
to 5 cs —30° in six bands (g p , r p , i p , z p , y p and w p ; 
iTonrv et al.ll2012f) to a limiting magnitude between 22.0 
mag to 24.0 mag depending on the filter, re-visiting a 
field e very 4 to 5 days. Fi nally, the Dark Energy Survey 
IDES: [Abbott et aLfl2005ft is observing part of the south¬ 
ern sky in four bands ( griz ) to a limiting magnitude of 
r ~ 24 mag, performing a roughly 30 deg 2 time domain 
survey for supernovae re-visiting a field every few days 
(4 to 5 days). 

One of the advantages of the QUEST-La Silla AGN 
Survey over these surveys is the very intense monitoring, 
observing the survey fields every possible night (although 
with large observing gaps from 2010 to 2012 - see more 
details in Section [3] below). We obtain between 2 to 5 
observations per night to remove spurious variability due 
to artifacts, to potentially study intra-night AGN vari¬ 
ability, and to produce stacked images to reach deeper 
magnitudes. Individual images reach a limiting magni¬ 
tude between r ~ 20.5 mag and r ~ 21.5 mag for a 
exposure time of 60 seconds or 180 seconds, respectively. 
One of the main characteristics of the QUEST-La Silla 
AGN Survey is that is focused on deeply observed extra- 
galactic fields with multi-wavelength coverage, and hav¬ 
ing nearly simultaneous observations in the near-infrared 
(near-IR) performed by the Visible and Infrared Survey 
Telescope for Astronomy (VISTA) surveys (see Section 
[3] below). 

X-ray variability studies indicate that the X-ray short- 
timescale normalization of the power spectrum den- 
sity (PSD) is correlated with the black hole mass (see 
McHardv 120131; IKellv et al .1120131 IGierliriski et al .1 [20081 
McHardy 198<f), and also that the black hole mass 
is correlated with the break frequency (char acteristic 
timescale) of the PSD dMcHardv et al.l 1200611 . These 
results have lead to the suggestion that AGNs are 
scaled versions of Ga lactic X-ray binaries with super- 
massive black holes ('Mc Hardv et all I2006T) . and that 
the X-ray characteristic timescale and the black hole 
mass are tightly related. Thus AGN variability could 
be a robust black hole mass estimator for a signifi¬ 
cant number of AGNs if these relations can be ex¬ 
tended to the optical/ultraviolet (UV) wavelengths. For 
a long time a significant anticorrelation between AGN 
variability amplitude and luminosity has been recog¬ 
nized in the opt ic al/UV regime (lUomoto et all 119761 
iHook et all 11994 iTrevese et al.1 11994 iCristiani et al l 


19971 IVanden Berk et all 12004 IWilhite et all 120081 : 
MacLeod et al.l 120101 IMeusineer et all 1201311 . Be¬ 
sides, some correlation between characteristic variabil¬ 
ity ti mescales and black hole mass have already been 
found (IKellv et al1l2009t [Collier fc Petersonl 1200 ill . Ad¬ 
ditionally, correlations between the optical variability 
amplitude and th e black hole mass have been described 
(Wol d et all 120071: IWilhite et al.l 120081 IKellv et al.l 120091 

IMacLeod et all 1201(1 IMeusineer et al l 1201311 . However, 
these results seem to point towards a more fundamental 
inverse relation between the AGN variability amplitude, 
at timescales longer than 1 year, and the Eddington ratio 
(or accretion rate) (Wilhite et all 120081: lAi et a. 1.1 [20101: 
MacLeod et al.l I201C : IZuo et al.1 l 2012 t IMeusineer et al l 
2013H . Further results indicate that AGN variability 
properties (e.g., the structure function) seem to change 
with other pro perties such as X-rav luminosity, and radio 
loudn ess (see IVanden Berk et ahl 12004 IMacLeod et al.1 
| 2010 L and references therein). All these results natu¬ 
rally lead to the conclusion that AGN variability is in¬ 
trinsically related to the physical parameters that govern 
the accretion in AGNs, and therefore, variability studies 
using large samples are crucial to improve our under¬ 
standing of the accretion process. 

One of the reasons why these relations are not com¬ 
monly used to estimate physical parameters is that until 
recently the sample of objects with good optical/UV/X- 
ray monitoring, having the required cadence of obser¬ 
vations to probe the necessary timescales ( up to years) 
was s mall. For exam p le, th e studies of iWold et al.1 
(1200711 and IKellv et al.1 (f200l were based on samples 
of only ~ 100 AGNs. Larger variability studies, on 


the other hand, hav 
identified AGNs on 

e been based on spectroscor 

hcally 

the Stripe-82 (IWilhite et al.l 

2008; 

MacLeod et al.l 20101 

Ai et al. 2010 : ; Schmidt et al. 

2010 ; 

Butler & Bloom 2011 

), which have ~ 60 photometric ob- 


servations over a deca de or on ensem ble AGN studies 
with few observations (IVanden Berk et~aLll2004l ). Since 
early AGN monitoring was mainly biased towards bright, 
color selected, or highly variable sources the calibration 
of the aforementioned relations is usually not represen¬ 
tative of the whole AGN population. It would also be 
highly desirable to compare black hole mass estimations 
obtained using different methods to assess systematic un¬ 
certainties (e.g., traditional methods vs. optical/X-ray 
variability black hole mass estimations). We expect to 
address most of these issues in the future using data col¬ 
lected as part of our survey. 

Since most of the AGN bolometric luminosity is emit¬ 
ted in the UV/optical part of the spectrum and their 
spectral energy distribution (SED), this is the region 
of choice to carry out an AGN survey. However, the 
UV/optical region can be subject to strong obscuration. 
On the other hand, the X-rays - particularly hard X-rays- 
can pass through the obscuring material, thus deep X-ray 
surveys provide a better census of the AGN population, 
particularly at low luminosities. X-ray observations are 
expensive and therefore surveys at these energies are usu¬ 
ally shallow or cover small areas. Clearly, to increase the 
synergy of any AGN variability survey it is necessary to 
carry out intense optical monitoring in fields with ex¬ 
tensive X-rays observations. This point is addressed in 
our survey by means of monitoring fields with extensive 
multi-wavelength coverage, ranging from the X-rays to 
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the mid-IR, and spectroscopic follow-up. 

Traditionally the optical selection of AGNs has made 
use of the fact that AGNs sh ow an UV excess in thei r 
SED compared to stars (see iSchmidt & Green! 11983 1. 
The UV-excess technique and the more recen t selectio n 
meth o ds ba sed on optical colors (see IRichards et al.l 
12002 L l2009h are very efficient in finding AGNs in 
the regions of color space where the AGN density is 
higher then the density of stars or galaxies. However, 
these selection methods based on optical colors miss 
a significant fraction of AGNs with peculiar colors 
(e.g., red QSOs) or QSOs located at a redshift range 
(2.5 < z < 3.0) wher e their opt i cal colors are similar 
to th ose of stars (see iFan 1 119991 : IRichards et al.l I2002L 
l2009f) . On the other hand, the fact that AGNs are 
highly variable makes their selection by means of 
variability a very promising technique to find them 
regardless of their colors. This selection method has 
been successfully used to identify a large number of new 
QSO c andidates dKozlowski et aD 120101: ISchmidt et al.1 


210101: iMacLeod et al. 20111: Butler fe Bloo 


Palana ue-Delabrouille et al. 12011: Kim et al 


am 

al. 


2011 


2012 


Graham^et^TL 2014 ^ _ Butler fc Blooml (|2011h and 


Palana ue-Delabrouille et al.1 (|201lfi have shown that 


selecting QSOs by their variability can increase consid¬ 
erably the number of QSO candidates in the redshift 
range where the colors of stars is similar to those of 
AGNs. 

Since 2010 we have been carrying out an AGN vari¬ 
ability survey using the wide-held QUEST camera lo¬ 
cated on the 1 m ESO-Schmidt telescope at La Silla 
Observatory. The telescope was fully robotized by 
Yale University, and the QUEST camera was moved 
from the 48” Palomar-Schimdt telescope to La Silla in 
2009. The aims of our survey are: 1) to test and im¬ 
prove variability selection methods of AGNs, and find 


niques (see Schmidt et al. 12010: 

Butler & Bloom 20111: 

iPalanaue-Delabrouillc et al. 2011 

); 2) to obtain a large 


number of well-sampled light curves, covering timescales 
ranging from days to years; 3) to study the link between 
the variability properties (e.g., characteristic timescales 
and amplitude of variation and other parametric vari¬ 
ability characterizations) with physical parameters of the 
system (e.g., black-hole mass, luminosity, and Eddington 
ratio). 

In this paper we present the technical description of 
our survey, and we study the relation of variability with 
multi-wavelength properties of X-ray selected AGNs in 
the COSMOS held. In Section^ we summarize the char¬ 
acteristics of the ESO-Schmidt telescope and the QUEST 
camera, and describe the operation of the robotic tele¬ 
scope. In Section [3l we present an overview of our ob¬ 
servation helds, and we discuss the reasons which make 
them special for carrying out an AGN variability survey. 
In Section [H we present the data reduction steps. In 
Section O we describe the steps to obtain the astromet¬ 
ric solution for the QUEST frames, and we present the 
quality of our astrometry. In Section E3 we dehne the 
Q-band photometric system, and describe in detail the 
steps to obtain well-calibrated light curves, including a 
nonlinearity correction to the photometry. In Section [Gj 
we also demonstrate the quality of our photometry, com¬ 
pare the quality of aperture and point spread function 


(PSF) photometry, and present examples of our light 
curves. In Section [7] we study the relation of optical 
variability with multi-wavelength properties of a sample 
of XMM-COSMOS X-ray selected AGNs. In Section [8] 
we summarize our work and present our conclusions. A 
standard A cold dark matter cosmology with Hq = 70 
km s -1 Mpc -1 , Dm = 0.27, and Ha = 0.73 is assumed 
throughout the paper. 

2. TELESCOPE AND CAMERA DESCRIPTION 

Schmidt telescopes are the instrument of choice for 
large area sky surveys because of their large field of view. 
The QUEST Camera was designed to operate at the 48” 
Samuel Oschin Schmidt Telescop e at the Palomar Ob¬ 
servatory Csee lBaltav et al.ll2007l) . Since 2009 the cam¬ 
era has been located in the 1-m Schmidt telescope of the 
European Southern Observatory (ESO) at La Silla, Chile 
(see TableQ]). This telescope is one of the largest Schmidt 
configurations in the southern hemisphere, situated in a 
dry site with dark skies and good seeing. Having a nearly 
identical optical configuration to the Palomar Schmidt, 
the QUEST camera was installe d at La Silla without an y 
changes to its front-end optics (Ra binowitz et al~ll2012fl . 
The surv ey uses the Q-band filter described below (Sec¬ 
tion 16.111 . 

Telescope pointing and camera exp osure are coo rdi¬ 
nated by a master scheduling program (Rab inowitz et all 
2012). A remote operator of one of the larger telescopes 
at the site (the ESO 3.6m) decides when conditions are 
appropriate for opening the telescope, and sends a re¬ 
mote command each night to enable the control software 
to open the dome. The control software automatically 
closes the dome whenever another nearby telescope (the 
2 .2m) is closed, when the sun rises, or when the remote 
operator sends a command to close. The remote oper¬ 
ator can monitor and control the state of the Schmidt 
telescope via a web-based interface. 

The CCD camera is located at the prime focus of the 
telescope about 3 m from the primary mirror. Its prop¬ 
erties are summarized in Table [lj The camera consists 
of 112 CCDs arranged in four rows or “fingers” of 28 
CCDs each as shown in Figure [U and covers 4.6° x 3.6° 
(north-south by east-west) on the sky. The fingers are 
flagged A, B, C, and D and the columns of CCDs from 
1 to 28. The gaps between the active areas of the CCDs 
in a row are typically 1 mm, and the gaps between the 
active areas in adjacent “fingers” are typically 22.8 mm 
corresponding to 25!8. Each CCD has 600x2400 pixels of 
13 /im x 13 /jm in size, and a pixel scale of 0'/882 pixel -1 . 
Therefore, to obtain a full coverage of the 4.6° x 3.6° field- 
of-view it is necessary to displace the telescope by half 
degree in right ascension to obtain an image-pair made 
of two tiles (see Figure [2J. 

Several CCDs have areas of high dark current owing to 
defective, electron-emitting pixels in the CCD substruc¬ 
ture or to a bad readout amplifier. By subtracting dark 
calibration images, the unfavorable influence of these de¬ 
fects on source detection is largely eliminated. However, 
about 16% of the CCDs are useless because they are per¬ 
manently off, randomly turn on and off, or have large 
defective areas which make it impossible to obtain an 
acceptable PSF, all of which hamper the astrometric so¬ 
lution due to the low number of stars detected and fake 
detections. The effective sky area covered by the func- 
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TABLE 1 

Properties of the ESO 1-m Schmidt 
telescope. 


Aperture Diameter 

1 m 

Focal Length 

3.05 m 

/-ratio 

3.05 

Plate scale 

14.74 /rm/ 

Latitude of Observatory 

-29° : 15' 

Longitude of Observatory 

o 

o 

Elevation 

2,375 m 


TABLE 2 

Properties of the QUEST Camera. 


Number ol CCDs 

112 

CCD Pixel Size 

13 /im x 13 fim 

Number of Pixels per CCD 

600 x 2, 400 

Pixel Size on Sky 

07882 x 07882 

Array Size, CCDs 

4 x 28 

Array Size, Pixels 

9,600 x 16,800 

Array Size, cm 

19.3 cm X 25.0 cm 

Array Size on Sky 

3.6° x 4.6° 

Sensitive Area 

9.6 square degrees 

Total Pixels 

161 x 10 6 



Fig. 1.— QUEST camera array, in black are shown death CCDs. 
The camera consist in four “fingers” of 28 CCDs. The fingers are 
flagged as A, B, C, and D, and the columns from 1 to 28. The 
gaps between adjacent fingers are 0.43°. Each CCD has 2400 X 600 
pixels of 0.59° X 0.15°. 

tioning CCDs is ~ 7.5 deg * 2 3 * * . 

3. SURVEY FIELDS 

We carried out the AGN variability survey observ¬ 
ing the COSMOS, ECDF-S, ELAIS-S1, XMM-LSS and 
Stripe-82 fields. These fields were chosen due to the 
wealth of ancillary data available for them, includ¬ 
ing XMM-Newton , Chandra , GALEX 1 HST , Herschel , 
Spitzer, and ground-based photometry and spectroscopy. 

Additionally, the COSMOS, ECDF-S, ELAIS-S1, and 
XMM-LSS fields have been repeatedly observed in the 
near-IR since 200 9, as part of the VIST A public sur- 
veys UltraVISTA (|McCracken et al.ll2012t l and VIDEO 

(j, Tar vis et all 1201,111 . The aims of these VISTA surveys 
are to study the evolution of galaxies out to z ~ 4, and 
achieve a comprehensive view of AGNs and the most 
massive galaxies up to the epoch of re-ionization. Ultra¬ 
VISTA observations were carried out over the same pe¬ 


riod as our survey, and an article presenti ng the stu dy of 
the A GN near-IR light curves is underway (1 Sanchez et al.1 

HSU). 

The equatorial Strip e-82 is part icularly valuable for 
variability studies (see iSesar et all 1200 7). This region 
was repeatedly observed as part of the Sloan Digital Sky 
Survey SDSS since 1998, this will allow us to extend light 
curves of some of Stripe-82 AGNs combining SDSS and 
QUEST photometry to more than 15 years. 

In Figures [2] and [3] we show the layout used to cover the 
COSMOS, ECDF-S, ELAIS-S1 and XMM-LSS fields. 
For the COSMOS, ECDF-S, ELAIS-S1 and Stripe-82 we 
used two tiles to cover the entire regions, while to cover 
XMM-LSS field we used four tiles (see Figure [3]). 

A brief summary of the observations is as follows. In 
2010 we obtained just a few observations for some fields 
(~ 20 nights). During 2011-2012, due to problems with 
the dome wheels of the telescope, we observed between 
100 to 150 nights per field (twice per night). From 2013 
onwards the observations have been performed more reg¬ 
ularly using an exposure time of 180 seconds (twice per 
night) reaching a magnitude limit of Q ~ 20.5 (r ~ 21), 
and observing each of our fields in more than 100 nights 
per year. Between March of 2011 and the end of Decem¬ 
ber 2014 we observed COSMOS on ~ 370 nights, ECDF- 
S on ~ 290 nights, ELAIS-S1 on ~ 500 nights, Stripe-82 
on ~ 450 nights and XMM-LSS on ~ 360 nights. In the 
case of XMM-LSS we usually observed 2 tiles per night, 
either tile 1 and tile 2 or tile 3 and tile 4. Observations of 
ECDF-S and XMM-LSS until 2012 have been reduced, 
while ELAIS-S1 and Stripe-82 are reduced until 2014. 
We expect to continue collecting data until mid 2016. 

From Section [5] onwards, the analysis presented corre¬ 
sponds to data obtained from 2010 to 2012 in the COS¬ 
MOS field. 

4. DATA REDUCTION 

For a typical night we take darks of 10, 60, and 180 
seconds, as well as morning and evening twilight flats. 
The exposure times for our science images are either 60 
or 180 seconds, with two image-pairs obtained per field 
per night. We reduce the QUEST data using our own 
custom IRAF0 scripts that carry out the reduction steps 
described below. 

Darks are combined to obtain a master dark for each 
detector and each exposure time. Then we subtract the 
master dark of the appropriate exposure time from the 
science images. The pixel-to-pixel variations of the sci¬ 
ence frames are corrected dividing them by the master 
flats obtained as we describe below. Finally, bad pixels 
in the science frames are interpolated using the IRAF 
task FIXPIX, which requires a bad pixel mask as in¬ 
put. We constructed bad pixel masks for each detector. 
To construct the bad pixel mask we use the IRAF task 
CCDMASK, taking as input the ratio between a high 
count twilight flat and a low count twilight flat. 

Given the large size of the camera and the large field of 
view of the telescope, to obtain dome flats using a uni¬ 
formly illuminated panel was unfeasible. Furthermore, 

4 IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 2.— Observing layout for the COSMOS (left) and ECDF-S (right) fields. In blue it is shown COSMOS1 and ECDF-S 1 (first tile) 
and in grey COSMOS2 and ECDF-S 2 (second tile). The second tiles are displaced by half degree from the first tiles. In red we show the 
area repeatedly observed by ultraVISTA and VIDEO in the COSMOS and ECDF-S fields, respectively. We flag the position of the four 
(A, B, C, and D) fingers for the first tiles. 




Fig. 3.— Observing layout for the ELAIS-S1 (right) and the XMM-LSS (left) fields. In blue are shown ELAIS-S1 1 and XMM—LSS1 1 
(first tile), in grey ELAIS-S1 2 and XMM—LSS1 2 (second tile), in magenta XMM—LSS2 1 (third tile), and in cyan XMM—LSS2 2 (fourth 
tile). Each tile is displaced by half degree from the previous. In red we show the area repeteadly observed by VIDEO, and we flag the 
position of the four (A, B, C, and D) finguers for the first tile. 


the automatic operation of the telescope complicates a 
careful acquisition of dome flats. Instead twilight flats 
were preferred because they can be easily acquired auto¬ 
matically. Typically twilight flats are taken at the begin¬ 
ning (evening flats), and the end (morning flats) of the 
observing night. On many occasions flats were not useful 
because of moon illumination gradients, a closed dome, 
or cloudy conditions. 

To avoid these problems and to obtain a good pixel- 
to-pixel variation correction we selected, from dark sub¬ 


tracted and trimmed twilight flats, the flats with mean 
counts above the average obtained over approximately 
two weeks of observations. Then we median combine 
these twilight flats and normalize them to obtain master 
flats for each detector. Typically we combined more than 
50 twilight flats per detector. 

5. ASTROMETRY 
5.1. Astrometric Solution 
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The frame headers lack basic astrometric information. 
To solve for this, we have created basic astrometric head¬ 
ers which contain a crude approximation to the true 
World Coordinate System information (WCS) of the 
frame. We introduce the estimated header keys (i.e., 
an estimated WCS) into the image headers using Miss- 
FITffl Then we proceed to obtain a first g ood astro- 
metric solution using astrometry. ne10 (lLang et al.ll2010l) . 
This software provides a fast and robust method to cali¬ 
brate astronomical images, and is able to produce a blind 
astrometric calibration of the QUEST frames. The out¬ 
put is an image with WCS information which can be 
further refined. 

The solution produced by astrometry.net is refined us¬ 
ing the Software for Calibrat ing AstroMetry and Pho¬ 
tometry (SCAMP, |Bertinj[2006|). SCAMP is a program 
that computes precise astrometric projection parameters 
from source lists obtained directly from FITS images. 
The input lists for S CAMP are in SExtractor binary for¬ 
mat (“FITS_LDAC”, iBertin fc Amoutsll 199 fill , and must 
contain the centroid coordinates, centroid errors, astro- 
metric distortion factors, flux measurements, and flux er¬ 
rors. Furthermore SCAMP creates frame headers ready 
to be used in an image stacking process. Both astrome- 
try.n et and SC AMP used the reference catalog USNO-B 
dMonet et al.l 120031 ) ■ which has an astrometric accuracy 
of 0"2. 


5.2. Assessment of QUEST astrometry 

To test the precision of our astrometric solution we 
assessed the internal astrometry by computing the stan¬ 
dard deviation of the internal cross-matching of sources. 
To assess the absolute astrometric precision we cross- 
matched our source catalogs with SDSS sources, and then 
we computed residuals of the cross-matching. 

5.2.1. Internal QUEST Astrometric Consistency 

To study the astrometric consistency we defined two 
samples of point-like sources: the subsample called stan¬ 
dard stars and the subsample of variable objects. We 
used a Xlc statistics (where LC s tands for light curve , 
generated as detailed in Section RT771 see lSesar et al.ll2007f l 
to distinguish between standard stars and variable ob¬ 
jects, which is defined as follows: 


2 1 ( m i “ w ) 2 

XLC ~(n~l)h 


( 1 ) 


where n is the number of detections, mi is the magnitude, 
m is the mean magnitude, and eq is the photometric er¬ 
ror. The standard stars subsample consists of objects 
with n> 10, a standard deviation in its light curve {ctlc) 
smaller than the median error observed in the photom¬ 
etry (er median), &lc < 0.035 mag, and xlc < L5 - The 
other group called variables contains objects with n > 5, 
ctlc > 2 x <T med i on , a LC > 0.05 mag, and xlc > 2 -°- As 
an example we obtained roughly 19,800 standard stars 
in our COSMOS field, of which 18, 265 have SDSS coun¬ 
terparts with good photometry. Of the 2, 790 variables 
objects in our COSMOS field roughly 2, 500 have a SDSS 
counterparts with good photometry. 


5 http://astromatic.net/software/missfits 

6 http://astrometry.net/ 



Fig. 4. — Astrometric comparison between QUEST and SDSS 
positions for all point sources in our COSMOS field. The figure 
shows the contour level lines of the number objects as a function of 
ARA and ADec, where ARA corresponds to &SDSS ~ &QUEST , 
and ADec to Ssdss ~ & QUEST- The black dashed lines show the 
zero offset axis, and the black pentagon marks the position of mean 
offset. 

We computed the mean and the median dispersion in 
the Equatorial (J2000) position of the standard stars and 
obtained 0/08 and 0'/07, respectively. The standard devi¬ 
ation for this subsample was 0/04. For the variable sub¬ 
sample the mean and the median were 0/15 and O'/11, re¬ 
spectively. The standard deviation for the variables was 
0/10. The variable subsample has many bright variable 
stars, and relatively faint transients. When the bright 
variable stars brighten, the PSF of these bright stars 
is close to saturation, such that the centroid becomes 
more uncertain. Something similar happens with faint 
transients where the centroid is more uncertain due to 
background fluctuations. In addition, after inspection of 
the light curves we found a significant number of non¬ 
variable sources classified as variables due to matching 
with nearby artifacts, thus leading to an increase in the 
astrometric dispersion. Overall our internal astrometric 
precision is typically ~ O'/I. 

5.2.2. QUEST vs. SDSS Astrometric Consistency 

To assess our overall accuracy we cross-matched our 
QUEST subsample catalogs with SDSS stars. For the 
standard subsample we obtained a mean of O'!21 with a 
standard deviation of 0'/13, and the median was 0/18. 
In the case of the variable group we obtained a mean of 
0'/23 with a standard deviation of O'/17, and a median 
of 0/19. As before some of the variable sources may be 
contaminated by matching with artifacts thus increasing 
the dispersion in their positions. 

In Figure 0] we show an astrometric comparison be¬ 
tween QUEST positions and SDSS positions for all point 
sources in our COSMOS field. We found mean offsets 
of asDSS — oiquest = 0/08 with a dispersion of 0/25, 
and Ssdss — Squest = —0/12 with a dispersion of 0/19. 
We found that these offsets are variable across the field 
(see Figure ©, and possibly are related to small sys¬ 
tematic differences between SDSS and USNO-B, since 
both astrometr y.net an d S CAMP used as reference cata¬ 
log USNO-B (|Monet et alJ l20031 which has an accuracy 
of 0/2. Therefore, the differences in the positions be¬ 
tween QUEST and SDSS objects are within the expected 
astrometric accuracy of the reference catalogs. 
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Fig. 5. — Astrometric comparion between QUEST and SDSS 
positions for twelve randomly selected detectors. In black we show 
tile one and in grey tile two. The black dashed lines show the 
zero offset axis. Overall the mean offsets are within 0.2" which is 
consistent with the astrometric accuracy of USNO-B llMonet et alj 

120031 . 

6. PHOTOMETRY 

6 .1. Q-band Photometric System 

The QUEST-La Silla survey uses a broad filter covering 
from 4000 A to 7000 A called the Q-band. This band¬ 
pass was designed to avoid the fringing often present in 
t he im a ges taken as part of the Palomar-QUEST survey 
(iBaltav et al.l 120071) . In Figure [6] we show an estima¬ 
tion of the system effective response of the Q-band (con¬ 
tinuous black line). We computed the Q-band system 
response profile using the QUE ST filter (dashed black 
line, see iR.abinowitz et al.l 120120 . the mirror reflectivity 
(dashed grey line), the Cerro Tololo Inter-American Ob¬ 
servatory (CTIO) sky transmission at an airmass of 1.3 
(dot-dashed grey line), and the quantum efficien cy of the 
camera (dotted grey line, see IBaltav et al.1^007 ). We as¬ 
sume a flat throughput of 90% for the optics of the cam¬ 
era. For comparison we show the filter response curves of 
SDS80 at an airmass of 1.3. It can be observed that the 
Q-band system response used in La Silla-QUEST is simi¬ 
lar to a broad ( g+r)sDSS filter. Therefore, we decided to 
calibrate the Q-band performing differential pho tometry 
using reference stars from SDSS-DR7 (lAbazaiian et, al.l 
1200911 in the COSMOS field. 

To calibrate our photometry we created a catalog of 
Q magnitudes of SDSS point sources with good g and r 
photometry. To make this catalog we transformed g and 
r magnitudes to flux, and then we added these fluxes 
to obtain an equivalent to the flux in Q. This flux was 
transformed to Q magnitudes in the AB system. 


6 .2. PSF Photometry 

We carried out PSF photometry in the Q UEST fr ames 
using custom scripts that run DAOPHOT ( Stetsonlfl987D 
on each epoch. We decided to compute PSF photome¬ 
try, since for faint stellar sources and crowded fields PSF 
photometry is usually better than aperture photometry. 
However, PSF photometry is more expensive in comput¬ 
ing time. Besides, we find that occasionally QUEST 
frames can yield a bad PSF due to telescope tracking 
problems, bad weather, moonlight scattered in to the 
frames, or frames with too many cosmetic problems. 


7 http://classic.sdss.org/dr7/instruments/imager/filters/ 


Fig. 6. — System responses of the SDSS bandpasses (in colors) 
and the Q-band (continuous black) at an airmass of 1.3. We com¬ 
puted the Q-band system response profil e from the multiplicatio n 
of the QUEST filter (dashed black, see IRabinowitz et al.l 120121 1. 
twice the mirror reflectivity (dashed grey), the sky transmission at 
an airmass of 1.3 (dotted-d ashed grey), and t he camera quantum 
efficiency (dotted grey, see IBaltav et af]|2007 ). We assume a flat 
throughput of 90% for the optics of the camera. The Q—band filter 
covers a spectral region similar to a (g + r)sDSS filter. 

6.3. Aperture Photometry 

We al so carried out aperture photometry using SEx- 
tractor (iBertin fe Arnoutii 119961) . SExtractor is a soft¬ 
ware to compute aperture photometry, and other source 
parameters (e.g., stellarity, source shape, and photome¬ 
try quality) quickly. 

To obtain AGN aperture photometry with minimal 
contribution from the host galaxy we determined an op¬ 
timal aperture of ~ 7 pixels (~ 6"18). This aperture 
is 2.5 to 3.0 times wider than the typical seeing of the 
QUEST-La Silla images, which ranges between 2"0 and 
2" 5, and is relatively insensible to typical seeing varia¬ 
tions at La-Silla observatory site. 

6.4. Comparison between Aperture Photometry and 
PSF Photometry 

Overall we found that PSF and aperture photometry 
agree very well and the differences are within the un¬ 
certainties. The mean difference between aperture and 
PSF photometry is 0.02 mag, with a standard deviation 
of 0.03 mag. In Figure [7] we compare PSF and aperture 
photometry for variable stars in the QUEST-COSMOS 
field. 

After a random inspection of the light curves (see Sec- 
tion l().7l) we found that PSF photometry achieves slightly 
better results for faint point sources, whilst for the rest 
of the objects both techniques yield similar results. How¬ 
ever, currently our pipeline fails to obtain a PSF solution 
for a significant number of image frames. Hence, for the 
AGN variability analysis of Section [7] we use aperture 
photometry. In the future we might be able to improve 
our PSF photometry using ad hoc parameters for each 
detector. 


6.5. Linearity Correction 

To ensure that non-linearity is not a dominant effect 
in the calibration of our data, we looked at observations 
of 60 seconds and 180 seconds taken close in time and on 
clear nights. In the perfect case where non-linearities and 
saturation were unimportant we would expect to see flat 
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Fig. 7. — Examples of periodic variable stars in the QUEST-COSMOS field, the left panel star has a period of ~ 1.44 hrs, while the right 
panel object has a period of ~ 14.59 hrs. We present PSF and aperture photometry in the top and middle panels, respectively. In the 
bottom panel we show the residuals of the difference between PSF and aperture photometry. 





Fig. 8. — Residuals of the difference between 180 seconds instrumental magnitudes (qisos) and 60 seconds instrumental magnitudes (960s) 
as a function of qisos for the detectors of Figure [9] The 180 seconds and the 60 seconds observations were taken on the same clear night 
5 min apart. These examples show how each detector has its own non-linear behavior. We plot black horizontal dashed lines to guide 
the eye. 
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Fig. 9. — Residuals of the difference between the QUEST instrumental magnitudes ( q ) and the ( g -\-r)sDSS magnitudes for stars plotted 
as a function of q. The black points correspond to the data used to compute the photometric zero-point as function of q (dashed-black 
line). Grey triangles are point sources with q <17.5, or q >13.5, or are outliers to the fit and they were not used to compute the best linear 
fit. On the top of each panel we give the best fit parameters and the detector name. 
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residuals for the difference between 180 seconds instru¬ 
mental magnitudes (gisos) and the 60 seconds instrumen¬ 
tal magnitudes (qeos)- In Figure [8] we show the residuals 
of < 7 i 80 s — 960 s as a function of gisos- In the plots all 
detectors show a clear upward trend toward the left cor¬ 
responding to bright objects. After an inspection of the 
images we found that the vast majority of the objects 
were not saturated, and the most likely explanation for 
this behavior is non-linearity of the detector. 

Furthermore, the onset of the non-linearities depends 
on each detector. For example, in Figure [8] detector A20 
has no significant slope, whilst detector C03 shows a clear 
positive slope. A pronounced slope is indicative that non¬ 
linearity is present. 

To correct the non-linearities we inspected the residu¬ 
als of the difference between the q instrumental magni¬ 
tudes and the (g + r)sDSS magnitudes as a function of q. 
In Figure [9]we show some examples of 9 — (g + r)sDSS as 
function of the instrumental magnitude. For most of the 
examples a linear fit as function of 9 is significantly bet¬ 
ter representation of the residuals than a constant zero 
point. For this reason we decided to use a linear expres¬ 
sion of the following form to fit the residuals: 

9 - (g + r)sDSS = a + /3 x q. (2) 

To avoid distortions in the fit due to very bright and 
faint sources we only used objects in the interval 13.5 <9 
<17.5 (shown in black), and we sigma-clipped the fit to 
eliminate outliers (grey triangles). Figure |9] shows that 
this model produces good fits to the residuals, even to 
the data points that were not used to compute the linear 
model. 

We determine the reduced y 2 (y 2 ), the standard devi¬ 
ation (<t), and the number of stars used to compute the 
fit. Subsequently we use these parameters to clean our 
photometry from bad epochs when constructing the light 
curves. We discard photometry with y 2 >10, er>0.1, or 
photometry that was obtained with less than 20 stars in 
the calculation of the zero-point, (see Section 16. 71) . The 
final calibrated Q-band magnitude is computed as: 

Q = 9 x (1 - /3) - a. (3) 

In Figure [TO] (top-panel) we show the residuals of Q 
minus the ( g+r)sDSS magnitude for 18, 265 non-variable 
stars where the calibration of the Q magnitudes was done 
following the prescription described above. These stars 
were classified as standards following the same criteria 
of Section [5.2.11 namely n > 10, a standard deviation in 
its light curve smaller than the median error observed in 
the photometry, a (Jlc < 0.035 mag, and Xlc < 1-5- As 
shown in the top-panel of Figure [TO] we obtained a good 
correspondence between {g + r)sDSS and Q magnitudes, 
with a mean residual of 0.004 ± 0.048 magnitudes. 

If instead we had used a constant zero point to cali¬ 
brate the observations, where we used cr-clipping to clean 
for outliers, we would have obtained 16, 977 standard 
stars with the ( g + r)sDSS ~ Q residuals shown in the 
bottom-panel of Figure [KJ] The correspondence between 
(9 + i~)sdss and Q is good with a mean residual of 
0.001 ±0.052. 

Using a linear fit instead of the usual constant zero- 
point reduces the dispersion in the light curves, thus in- 



Q 


Fig. 10.— Contour level lines of the number of stars as a func¬ 
tion of the residual magnitude, and the Q-band magnitudes for 
more than 15, 000 stars. The residual is the difference between the 
(g + r)sDSS magnitude and the mean Q-band magnitude of the 
stars. These stars were flagged as standards based on our 2011 
observations. In panel (a) the Q-band calibr ation was done using 
our linearity correction technique (see Section l6.5ll . whilst in panel 
(b) the Q-band calibration was done using a constant zero point. 

creasing by ~ 7.5% the number of non-variables or stan¬ 
dard stars. Although the means and the standard devia¬ 
tions of the (g + r)sDSS — Q distributions in the case of 
a linear and constant zero-point are consistent with each 
other, a Kolmogorov-Smirnov (K-S) test yields a p-value 
equal to 0.01. This means that the null hypothesis that 
both distributions are drawn from the same parent pop¬ 
ulation can be rejected (see Figure [FT]). Additionally, 
a K-S test on the Q magnitude distributions yields a p- 
value lower than 0.01, so the null hypothesis that both Q 
magnitude distributions are drawn from the same parent 
population can be rejected. 

We found that when computing the zero-point using 
a linear model with Eqs. [2] and [3] the number of non¬ 
variables or standard stars increases significantly in the 
range 15 <Q <19 (Figure [Till . At the same time the 
linear model keeps the mean and the standard deviation 
of the (g + r)sDSS — Q residual consistent with the case 
of a constant zero-point, and reduces the skewness in 
(9 + r)sDSS — Q as function of Q (see Figure [TUI). 

In summary, we found that the use of a linear model 
to compute the zero-point as a function of 9 achieves a 
better calibration of our data, and in the cases where 
non-linearity is unimportant (i.e., small /3) the photom¬ 
etry remains consistent with the case of a constant zero- 
point. 

6 . 6 . Color-term 

The agreement between between the (9 + r)sDSS and 
the Q-band is overall very good. However, we decided to 
explore the possibility of adding a color-term to enhance 
our calibration. 

With an estimation of the Q-band system response at 
hand, we looked for a relation between the residuals of 
Q — (9 ± r)sDSS as a function of the (9 — r)sDSS color. 
We began by multiplying t h e spectroph otometric stan¬ 
dard stars of IHamuv et all (11994 1199211 with the sys¬ 
tem response curves of Figure [ 6 ] Then we computed 
Q - (9 + r)sDSS , and the (9 - r) S DSS sy nthetic co lors. 
The results of the synthetic photometry of lHamuv et al.1 
(1994, 1992) standards are shown as black circles in Fig¬ 
ure [I2j We modeled the Q — {g + r)sDSS residuals as a 
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Fig. 11. — Left panel: normalized cumulative distributions of 
Q — (g-\-r)sDSS residuals for non-variable point sources for the case 
of a constant zero-point (grey), and for a linear model calibration 
(black). Right-panel: normalized cumulative distributions of Q 
magnitude for non-variable point sources for the case of a constant 
zero-point (grey), and for a linear model calibration (black). The 
p-value and the Kolmogorov-Smirnov D parameter are indicated 
on the left of each panel. 
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Fig. 12 . — Q — (g-\-r)sDSS residuals of the synthetic photometry 
versus the (g — r) spss syntheti c color s of the spectrophotometric 
standard stars of IHamuv et all < 119941 , 119921 ) . The black line is 
the best fit to the stars. Additionally, for comparison we show 
the Q — (g + r)sDSS residuals of the synthetic photometry versus 
(g — r ) SPSS synthetic colors of t he composite quasar spectra from 
SDSS dVanden Berk et al.l 120011 ) , in addition to A GNs, galaxies 
and s tarburst galaxies templates from SWIRE (see PoHol la et all 

|200W) . 

function of the ( g — r)sDSS colors using a linear relation 
(dashed line in Figure fl2l) of the form: 


Q ~ ( g+r)sDSS = —0.359(±0.001) -0.0686(±0.0014)x( 5 

(4) 

To estimate the uncertainty in the parameters of the fit, 
we selected randomly 20 of the 29 stars and then we re¬ 
fitted a linear model. We repeated this random selection 
50 times, then we estimated the standard deviation of 
the adjusted parameters. 

We also obtained synthetic photometry of extragalac- 
tic sources to investigate if they follow a relation similar 
to that of the stars. In Figure[[2]we show synthetic pho¬ 
tometry obtained from SWIRE templates of AGNs, ellip¬ 
tical galaxi es, sp iral galax ies, an d starburst galaxies (tri¬ 
angles; see lPolletta et al.l I2007T) . We additionally show 
synthetic photometry obtained from a SDSS compos¬ 


ite qu asar spectrum (blue pentagon; IVanden Berk et ahl 
1200111 . Overall, the template of extragalactic sources fol¬ 
lows a linear relation similar to that of the stars. The 
outliers are AGNs, and they depart from the main rela¬ 
tion mostly due to the presence of strong emission lines. 

We used Eq. 0] to transform the photometry from the 
(g + r)sDSS photometric system to the Q-band. Then 
we tested our transformation performing differential pho¬ 
tometry on several nights of 2011. We found that defin¬ 
ing the Q-band as in Eq. d] does not produce better 
results than defining the Q-band as equivalent to the 
(g + r)sDSS- Indeed using Eq. []]to compute the Q-band 
magnitude for the reference stars and then to compute 
the zero points, increases the dispersion in the zero point 
calculation and overall produces worse fits. 

As a second check on the effects of adding a color- 
term, we performed differential photometry for reference 
stars on four clear nights of 2011. We then computed the 
observed color-term fitting a straight line to the differ¬ 
ence between the instrumental magnitude ( q ) minus the 
reference ( g + t)sdss magnitudes as a function of the 
{g-r) S DSS colors. In Figure H31 (left panel) we show the 
distribution of the observed color-terms obtained. The 
mean of the distribution is 0.004, the median 0.007, and 
the standard deviation 0.028. These values are consistent 
with no color-term. However, the large dispersion in the 
color-term distribution may be explained by a color-term 
variation between CCDs. Small variations from night-to¬ 
night are also expected. 

Formally the difference between the computed color- 
term using the system responses and the observed color- 
term is 2.59 a. The discrepancy is most likely due to 
differences between the estimated and the real Q-band 
system responses. For example, we made a rough as¬ 
sumption of the optics transmission of the camera, and 
almost certainly there are differences in the quantum ef¬ 
ficiency from detector-to-detector as function of wave¬ 
length. 

In conclusion we find that observations indicate that 
the color-term is close to zero, and therefore the Q-band 
is well described by the ( g + t)sdss system. As result 
we did not apply any color-term to our transformation. 
Nevertheless, there is some uncertainty associated to this 
assertion. In order to correctly take into account this 
uncertainty we added in quadrature a color-term uncer¬ 
tainty equal to 0.03 x (g — r)sDSS- 

In Figure [13] (right panel) we show the resultant dis¬ 
tribution of the observed color-terms after adding this 
color-term uncertainty. The mean of the distribution was 
0 .00002, the median was 0.002, and the standard devia- 
—r)^flg^ibtained was 0.029. Overall this additional term did 
not increase the dispersion in the zero point calculation 
and produced consistent fits. Therefore we assume that 
it is a robust measure of our uncertainty on the filter 
transmission discrepancies. 

6.7. Light Curve Generation 

For each epoch and each detector we generated a cata¬ 
log which contains the Equatorial position of the object 
in degrees (i.e. right ascension and declination), the cal¬ 
ibrated magnitude, the magnitude error, the number of 
stars used to obtain the photometric zero-point, the stan¬ 
dard deviation of the zero-point, and the reduced y 2 of 
the zero-point adjustment (see Section [63]). Light curves 
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Fig. 13. — Distribution of the color-terms computed from ob¬ 
servations of four nights made on 2011. (Left panel) color-term 
calculated from a Q-band catalog of stars computed assuming a 
equivalence between the Q-band and the (g + r)s£),ss-band (i.e. 
zero color-term) but not taking into account the uncertainty in the 
color-term correction. (Right panel) color-term calculated from a 
Q-band catalog stars computed assuming a equivalence between 
the Q-band and the (g-\-r)sDSS~ band but taking into account the 
uncertainty in the color-term correction (see text). The mean, the 
median and the standard deviation of the distributions are given 
on the top-left of each panel. 


TABLE 3 

Magnitude Difference of the Overlapping Stars. 


Range 

N Stars 

<A Q> 

(7 

AQ Median 

Q <14.0 

40 

-0.044 

0.065 

-0.036 

14.0 <Q <15.5 

316 

-0.020 

0.056 

-0.016 

15.5 <Q <17.0 

1501 

-0.019 

0.068 

-0.015 

17.0 <Q <18.5 

2844 

-0.015 

0.067 

-0.013 

Q >18.5 

1346 

-0.014 

0.073 

-0.014 

15.5 <Q <18.5 

4345 

-0.016 

0.067 

-0.014 

Total 

6047 

-0.016 

0.068 

-0.014 


are generated for objects in the matching catalogs of dif¬ 
ferent epochs, but always for the same detector, and the 
same tile. This means that we generate light curves for 
each tile/detector separately. To match sources we used 
a radius of 1", which is roughly equivalent to 1 pixel. 
We avoid contamination from bad nights, which typi¬ 
cally have few stars in the zero-point computation and 
produce a very uncertain photometric zero-point. To this 
effect, we select only those nights that use more than 20 
stars to compute the zero-points, where xt < 10, and the 
standard deviation of the zero-point adjustment is lower 
than 0.10 mag. 

6 .8. Systematic Error in the Photometric Zero-Point 

We assessed systematic differences in our photometry 
by comparing the magnitudes of stars in the common 
area of two adjacent tiles offset by 0?5. We found that 
there are 6047 stars detected in two adjacent images of 
the COSMOS field. The mean difference in their mag¬ 
nitudes is 0.016 ± 0.068 magnitudes, with a median dif¬ 
ference of —0.014 magnitudes. In Figure [14] we summa¬ 
rize the mean magnitude differences for each overlapping 
area. We investigated the magnitude differences as a 
function of magnitude dividing the stars in magnitudes 
ranges, and we found that stars in the brighter range have 
higher systematic differences than fainter stars. Our re¬ 
sults are summarized in Table [I] 


Fig. 14.— Mean differences for the stars in common for each 
overlapping area. The corresponding overlapping area is indicated 
in the top-left of each panel. 

Based on our results we quote that our systematic 
uncertainty is 0.05 magnitudes, this value brings into 
agreement most of our observed/calibrated photometry 
in different tiles. We note that this systematic uncer¬ 
tainty is related to variations from detector-to-detector, 
and that additional variations are also introduced as a 
consequence of the use of different stars to compute the 
zero-points. For observations obtained using the same 
tile/detector it is not necessary to take in to account this 
uncertainty, since we are using the same detector and 
mostly the same stars to compute the zero point. 


7. OPTICAL VARIABILITY OF XMM-COSMOS X-RAY 
SELECTED AGNS 

In this section we study the optical variability proper¬ 
ties of the XMM-COSMOS X-ray selected AGNs. To 
do this we constructed optical light curves using the 
2010-2012 data for 287 QUEST detected sources out 
of the_1797 X-r ay point sources presented in Table 
2 of iBrusa et al.1 (|2010T) (the X-ray catalog hereafter). 
Most of the objects contained in the the X-ray cat¬ 
alog are faint optical sources with a mean r magni¬ 
tude of ~ 22.9 and a standard deviation of 1.9 mag 
in their r magnitude distribution. Overall, we detected 
~ 76% (254 objects) of the sources in the the X-ray 
catalog to a limiting magnitude of r ~ 21. The X-ray 
catalog contains X-ray, optical, infrared, spectroscopic, 
and photometric redshift information for the XMM- 
COSMOS sou r ces. To construct this multiwavelength 
data set IBrusa et al.1 (|2010f ) cross-correlate d X-ray posi¬ 
tions with the optical multiband catalog of lCapak et al.1 
(|2007f) . the Can ada-France-Hawaii Te lescope (CFHT) 
K-ba nd catalog d McCracken et al .ll2010D. th e IRAC cat¬ 
alog (iSanders et al. 2007f Ilbert, et al.1 200if). and the 24 
H m MIPS catalog ( Le Floc’h et al.ll 2009f). Additi onally, 
they used accurate C-Chandra (|Elvis et all 120091 posi¬ 
tions available for a subset of objects within an area of 
~ 1 deg 2 to control-check the optical/near-IR identifica¬ 
tions, and to assess the reliability of the proposed iden¬ 
tifications. 

Good-quality spectroscopic redshifts for the opti¬ 
cal counterparts are available from Magellan/IMACS 
and MMT observation campaigns (~ 530 objects; 

I Trump et al.1120071 f2009l). from the VIMO S/zCOSMOS 
)iects: lLi~ 


(~ 580 objects: ILillv et all 120071 l200l , or were al¬ 
ready present either in the SDSS catalog (~ 100 objects; 
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lAdelman-McCarthv et ai1l2006t iKauffma nn et al l 2003) 


or in the literature (~ 95 objects: iPrescott et alJ 2006 1. 

Using this large spectroscopic data set. iBrusa et alJ 
• 2010; divided the extragalactic sources with available 
spectra into three classes, on the basis of a combined 
X-ray and optical spectroscopic classification. 


• Broad-line AGN (BL AGN): all objects having 
at least one broad (FWHM >2000 km s -1 ) opti¬ 
cal emission line in the available spectrum; (421 
sources). 

• Non-broad-line AGN (NL AGN): all objects with 
unresolved, high-ionization emission lines, exhibit¬ 
ing line ratios indicative of AGN activity; if lines 
are not detected or the observed spectral range 
does not allow to construct line diagnostics, objects 
with rest-frame hard X-ray luminosity in excess of 
2 x 10 42 erg s _1 , typical of AGNs; (370 sources). 

• “Normal” galaxies (GALs for short): all sources 
with unresolved emission lines consistent with spec¬ 
tra of star-forming galaxies or galaxies showing 
only absorption lines, and with rest-frame hard X- 
ray luminosity lower than 2 x 10 42 erg s _1 , or un¬ 
detected in the hard band; (53 sources). 


I White dwarfs 

II Low-z QSOs 

III dM/WD pairs 
2.0 - IV rr Lyrae stars 

V Stellar locus stars 

VI High-z QSOs 



0 12 3 4 

u-g 


Fig. 15. — u — g vs. g — r color-color diagram for all the XSMM- 
COSMOS extragalactic sources detected in the QUEST-La Silla 
survey. Blue points correspond to BL AGNs, red triangles to NL 
objects, and green squares to GAL sources. Larger symbols are 
used to highlight the bes t sample, and white spot identify variable 
sources (see Section 17.21) . To guide the eye we show the division 
used in loesar et al.l (120071) as grey dashed lines, that it corresponds 
to regions where different subclasses of variables objects are dom¬ 
inant. The classes are listed on the top-left of the figure. The 
photometry correspond to model magnitudes obtained from the 
SDSS-DR7 database, using a cross-matching radius of 073 to match 
objects. The photometry was corrected using Galactic extinction 
values given by SPS S, that make use of the extinction maps of 
ISchlegel et al.l (| 19981 ). 


7.0.1. u — g color vs. g — r color 

In Figure fl5l we show the u — g vs. g — r color-color di¬ 
agram for all the XMM-COSMOS extragalactic sources 
detected in our survey construct ed using SDSS-DR7 pho¬ 
tometry ( A bazaiian et al.l 120091 ). In the figure blue cir¬ 
cles correspond to BL AGNs, red triangles to NL AGNs, 
and green squares to GALs. We use larger symbols to 
high light the best observed sample defined below in Sec¬ 
tion o an d a white spot to identify variable sources 
(see Section 1731) . To guide the reader’s eye we plot the 
division used in lSesar et al.l ()2007l ) as grey dashed lines, 
which corresponds to regions where different subclasses 
of variables objects are dominant. For example, in re¬ 
gions II and VI the dominant classes of variable objects 
are low-z and high-z QSOs, respectively. 

The photometry corresponds to model magnit udes ob¬ 
tained from the SDSS-DR7 database (lAbazaiian et al.l 
1111 ), using a cross-matching radius of 0'.'3 to match ob¬ 
jects. Using this cross-matching radius we find that 97% 
of BL AGNs, 88% of NL AGNs, and 83% of GALs have a 
SDSS counterpart. The photometry was corrected using 
Galactic extinction values given by S PSS, that make use 
of the extinction maps of ISchlegel et all (119981) . 

From the figure it is clear that in the u — g vs. g — r 
color-color space BL AGNs, NL AGNs, and GALs oc¬ 
cupy roughly different regions. While BL AGNs tend to 
be concentrated in region II (low-z QSOs), with a few 
objects in region VI (high-z QSOs), and a non-negligible 
fraction of objects in region III (dM/WD pairs); NL 
AGNs are concentrated in region III, with some objects 
in region V (Stellar locus stars), and very few objects 
in region II. On the other hand, GALs are concentrated 
in region V, with a non-negligible fraction of objects in 
region III. In particular, GAL objects seems to occupy 
the color-color region between region III and region V. 

As expected, BL AGN colors are dominated by AGN 
emission, while GAL colors are dominated by stellar 


emission. On the other hand, the optical colors of NL 
AGNs are consistent with obscured AGNs with some con¬ 
tribution from the host galaxy stellar emission. Since 
IBrusa et~all ( 2010 ) classification is partly based on opti¬ 
cal spectra this confirms that they classification is robust 
and consistent with the optical colors. 

7.1. QUEST-La Silla Variability Survey Completeness 
and Detection Limit 

Of the 287 XMM-COSMOS sources detected by 
QUEST, 233 sources have a corresponding r-band mag¬ 
nitude in the X-ray catalog. The remaining 44 sources are 
bright, probably saturated, and therefore not included in 
the optical multiband catalog. In Table [T] we summarize 
the detection completeness of our survey as function of 
r magnitude ranges for easy comparison with other sur¬ 
veys, we also breakdown our analysis for different classes 
of extragalactic objects. Table |T| is organized as follows: 
r magnitude range (Column 1); total number of sources 
detected in our survey in the corresponding magnitude 
range (Column 2); number of objects detected in our sur¬ 
vey over the total number of XMM-COSMOS objects in 
the corresponding magnitude range (Column 3); num¬ 
ber of BL AGNs detected in our survey over the total 
number of XMM-COSMOS BL AGNs in the same mag¬ 
nitude range (Column 4); number of BL AGNs detected 
in our survey over the total number of sources detected 
in our survey; in parentheses the number of total XMM- 
COSMOS BL AGNs over the total number of XMM- 
COSMOS sources in the magnitude range (Column 5); 
number of NL AGNs detected over the XMM-COSMOS 
NL AGNs in the corresponding magnitude range (Col¬ 
umn 6); number of NL AGNs detected over the total 
number of sources detected in the range, and in paren¬ 
theses the number of XMM-COSMOS NL AGNs over the 
total number of XMM-COSMOS sources in the range 
(Column 7); number of GALs detected over the total 
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number of XMM-COSMOS GALs in the range (Column 
8 ); number of GALs detected over the total number of 
sources detected in our survey, and in parentheses the 
number of total XMM-COSMOS GALs over the total 
number of XMM-COSMOS sources in the range (Col¬ 
umn 9). 

The survey is roughly 75% complete in the magnitude 
range between r ~ 17 mag and r ~ 20 mag, and slightly 
less complete (67%) up to r ~ 21 mag. In the range 
between r ~ 21 and r ~ 22 mag the survey declines to 
12 % completeness, and the survey drops to less than 1 % 
completeness at fainter magnitudes. If we break down 
the detection completeness fraction for different object 
classes (n^ e “ ss /ritot ss ; see Table [lj, we find that the de¬ 
tection completeness fraction, for all classes of objects, 
is consistent with roughly 75% up to r ~ 20 mag. In the 
r magnitude range between 20 mag and 21 mag the de¬ 
tection fraction of BL AGNs is 74%, whilst the detection 
fraction of NL AGNs and GALs is 60%. At fainter mag¬ 
nitudes the detection fraction of BL AGNs is significantly 
higher than that of NL AGNs and GALs. This could be 
due to strong UV rest-frame emission lines (Ly-a, Si IV, 
C IV, and C III) entering the blue side of the Q-band for 
sources at z >2.0. When we break down the number of 
objects detected of a certain class over the total number 
of objects detected in a magnitude range {n^ t ss /rideu see 
Table [l]), they are in very good agreement with the num¬ 
ber of XMM-COSMOS objects of the same class over the 
total number of XMM-COSMOS objects in the sample 
( n tot ss /ntot', see Table Q} in the magnitude range up to 
r ~ 21 . 

The analysis of the detection fractions indicates that 
we are detecting a large fraction of objects up to r ~ 21 . 
This is in agreement with the limiting magnitude ex¬ 
pected based on the exposure time used; for an exposure 
time of 60 seconds the expected limiting magnitude is 
Q ~ 19.5 which is equivalent to a r magnitude in the 
range between 20 mag and 21 mag. Similarly, for an ex¬ 
posure time of 180 seconds the expected limiting magni¬ 
tude is Q ~ 20.5 equivalent to a r magnitude in the range 
between 21 mag and 22 mag. During 2011 and 2012 most 
of our observations were performed using exposure times 
of 60 seconds, where the completeness below r ~ 21 mag 
is low. Currently (2013-2014) we are performing our ob¬ 
servations using exposure times of 180 seconds which will 
increase our completeness at fainter magnitudes. 

To illustrate our detection limit as function of redshift, 
in Figure fT6l we show the Q-band absolute magnitude of 
the extragalactic XMM-COSMOS objects as function of 
redshift. The Q-band absolute magnitudes are not K- 
corrected or corrected by dust extinction. As can be seen 
we are detecting objects up to redshift ~ 3. In total we 
detect 29 GALs all below redshift 0.4, 42 NL AGNs all 
below redshift 1.0, and 123 BL AGNs. In Figure [16] we 
additionally show the detection limit for an object with 
Q = 19.5 mag (continuous line), and the detection limit 
for an object with Q = 20.5 mag (dashed line). 

7.2. Variability of the XMM-COSMOS Objects 

To study the variability of the XMM-COSMOS ex¬ 
tragalactic sources we defined a flux limited sample (Q 
<19.5), where the time lapse between the first and the 
last observations was more than 600 days. We call this 
the best observed sample. The time span ensures that 



Fig. 16.— Q-band absolute magnitudes for XMM-COSMOS 
sources detected as part of our survey. Absolute magnitudes are 
not K-corrected or corrected by dust extinction. Blue points cor¬ 
respond to broad-line AGN, red triangles correspond to objects 
classified as non-broad-line AGN, and green squares are classified 
as “normal” galaxies or starburst galaxies. The continuous line is 
the detection limit for an object with Q = 19.5 magnitude, and 
the dashed line is the detection limit for an object with Q = 20.5 
magnitude. 

the light curves have a significant number of observa- 
tions. Additionally the time span longer than 600 days 
ensures that we will be able to study long-term variabil¬ 
ity as this corresponds to rest-frame timescales of at least 
a few months for objects at z > 0.6, and of > 1 year for 
objects at z < 0.6. Our sample of best observed light 
curves consists of 55 BL AGNs, 24 NL AGNs, and 23 
GALs. 

To classify an object as var iable or non-variable we 
used the variability index V (iMcLaughlin et al.1119961 
iPaolillo et, al.ll2004iLanzuisi et al.ll2014fl . defined as V = 
— log(l — P(x 2 )), where P{\ 2 ) is the the probability that 
a x 2 lower or equal to the observed value could occur 
by chance for an intrinsically non-variable source. We 
also defined P var = 1 — P{\ 2 ) as the p robability of an 
object of being variable. Following lLanzuisi et al.1 (|2014H 
we defined V = 1.3 as our threshold to define a variable 
object, since an object with V >1.3 has more than 95% 
probability to be variable ( P var >0.95). Based on the 
V index we find that 44 (80 ± 5%) of our best observed 
BL AGNs are variable, 5 (21 ± 8%) of our best observed 
NL AGNs are variable, and 5 (22 ± 9%) of our best ob¬ 
served GALs are variable. Since the division between 
variable and non-variable objects is somewhat arbitrary, 
we assumed a binomial distribution in the classification 
to estimate the errors. Therefore, based on the binomial 
distribution 3 BL AGNs, 2 NL AGNs, and 2 GALs may 
be classified as variable sources by chance. 

As a measure of the variability amplitude we used the 
excess variance defined as cr^ms = °lc ~~ a median > where 
<7lc is the standard deviation of the light curve, and 
a median is the median error of the light curve. In the 
case where a\ c <of n ed% an we used the median error as 
an upper limit on the excess variance. In Table Q] we 
present the distribution of a rms for our sample of best 
observed objects divided in three bins (oy ms < 0.05, 
0.05 < a rms < 0.10, and a rms > 0.10). We found 
that BL objects tend to have larger variability amplitude, 
with 77% of the objects showing a rms > 0.05, and 23% 
of the objects with a rms > 0.10. On the other hand, 33% 
of the NL AGNs and GALs objects have a rms > 0.05. 













14 


Cartier et al. 


TABLE 4 

QUEST-La Silla AGN Variability Survey Completeness in the COSMOS Field. 


Range 

n c£ei 

n det/ n tot 

n ! ’LMot 

n “ ,/ n det ( n tit/n tot) 

v~.nl /„nZ 
n det/ n tot 

n ™lt/ n det «it/nt ot ) 

qa / qa 
n d,J n tot 

n ®“,/ n det {'O-int/'O-tot) 

i 

2 

3 

4 

5 

6 

7 

8 

9 

17 <r <18 

9 

0.75 


0.00(0.00) 

1.00 

0.33(0.25) 

0.62 

0.56(0.67) 

18 <r <19 

26 

0.87 

0.80 

0.31(0.33) 

0.80 

0.15(0.17) 

0.88 

0.27(0.27) 

19 <r <20 

58 

0.74 

0.72 

0.45(0.46) 

0.73 

0.19(0.19) 

0.75 

0.16(0.15) 

20 <r <21 

107 

0.67 

0.74 

0.62(0.57) 

0.59 

0.16(0.18) 

0.60 

0.06(0.06) 

21 <r <22 

30 

0.12 

0.20 

0.67(o.40) 

0.06 

0.17(0.32) 

0.00 

0.00(0.04) 

22 <r <23 

3 

0.01< 

0.02 

0.67(0.33) 

0.01< 

0.33(0.33) 

0.00 

O.OO(O.OK) 


TABLE 5 

Excess Variance Distribution. 


Range 

BL 

NL 

GAL 

&rms <0.05 

0.23 (0.20) 

0.67 (0.42) 

0.67 (0.56) 

0.05 <cr rms <0.10 

0.54 (0.56) 

0.17 (0.50) 

0.33 (0.35) 

&rms >0.10 

0.23 (o.24) 

0.17 (0.08) 

0.00 ( 0 . 09 ) 


Note. — Fraction of BL, NL, and GAL objects in a given 
range of a Trns for the best observed sample. In parentheses 
we give the fraction of objects in a given range, but now 
including upper limits. 

Of these, 1 NL AGN (17%) have a rms > 0.10 while all 
GALs have a rms below this value. When we include up¬ 
per limits to the distributions the tendency remains the 
same with 80% of BL AGNs, 58% of NL AGNs, and 44% 
of GALs having oy ms > 0.05. 

7.2.1. Excess Variance vs. [3.6] — [8.0] Color 

The mid-IR colors have been propose d as a se¬ 
lection method t o identify AGNs (see lLacv et all 
gOO^ jSternel^l.l 120051 : iKozlowski fc Kochanekl 120091 : 
Assef et al.1 12013L and references therein). The method 
relies on the fact that the mid-IR emission of normal 
galaxies and AGNs have different spectral energy dis¬ 
tributions (SED). The composite spectra of the stel¬ 
lar population of normal galaxies produces a SED that 
peaks at approximately 1.6 /xm, while AGNs have a 
roughly power-law SED due to hot dust emission, that 
reprocess the UV/optical emission o f the accretio n disk, 
peaking somewh ere around 10 ^m dMor fe Netzeill2012l 
iLira et all 120131) . Following this idea lBrusa et al.l (1201011 
showed that the IRAC broadband [3.6] — [8.0] mid-IR 
color can be used to disentangle objects with infrared ris¬ 
ing SED (red power-law), e.g. AGNs, and sources with 
an inverted SED (blue power-law), e.g., normal galaxies 
at low-z. They used this method to select a sample of 
highly obscured luminous AGN candidates. 

In Figure [lT] we explore a tentative relation between 
the a rm s and the [3.6] — [8.0] mid-IR color for our sam¬ 
ple of best observed BL AGNs, NL AGNs, and GALs. 
In the Figure we show BL AGNs with z >1 using cyan 
circles, BL AGNs with 2 <1 using blue circles, red tri¬ 
angles correspond to NL AGNs, and green squares are 
GALs. The grey dashed line corresponds to the division 
between AGN dominated and host dominated sources of 
iBrusa et alJ ([20101 ) ([3.6] — [8.0] = 0.856), and upper lim¬ 
its are shown as grey downward arrows for BL and NL 
AGNs, and not shown for GALs for clarity. 

As can be seen in Figure[l7]BL AGNs show on average 
large a rms . Among them there is an apparent trend, 
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Fig. 17. — errms vs. [3.6] — [8.0] color for the sample of best 
observed XMM-COSMOS objects. In the figure cyan cirlces are BL 
AGNs with z >1, BL AGNs with z <1 are shown as blue circles, 
red triangles correspond to NL AGNs, and green squares are GALs. 
The grey dashed line corresponds to the division be t ween AGN 
dominated and host dominated sources of IBrusa et al.l ([20101 1. and 
upper limits are shown as grey arrows down for BL and NL AGNs, 
and not shown for GALs for clarity. 


although with large scatter, where AGNs with higher 
arms tend to have redder [3.6] — [8.0] colors. This trend 
is largely driven by BL AGNs, particularly those at z >1 
(cyan circles), where the expected contribution from the 
host galaxy to the Q-band is very low. Among the BL 
AGNs with z >1 the trend seems not to be related to 
redsliift or luminosity. Finally, the BL AGNs with z <1 
have bluer colors than the overall BL AGN sample which 
implies a higher contribution of host galaxy light to the 
photometry. 

A plausible interpretation for the trend shown in Fig¬ 
ure [EH is that AGN variability is diluted by host galaxy 
light in objects with bluer mid-IR colors (i.e., more host- 
dominated). This will be explored by our collaboration 
in the future using a larger sample obtained from pho¬ 
tometry of stacked images and with light curves with 
longer time span. 

7.2.2. Structure Function 

The structure function is a simple tool to quantify the 
variability of a source as a function of the time lapse 
between observations (r) when a significant number of 
observations are available. To study the variabilty be¬ 
haviour of our sample as a function of time we calculated 
the structure function as: 


SF(t) 
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i 

N 


N 

7=1 


m(ti + t)} 2 . 


(5) 
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where m(t; ) is the magnitude at time ti, m(ti + r) is the 
magnitude at time tj + r in the rest-frame, and N is the 
number of observations in the time span bin. To compute 
the structure function we average values in equal size bins 
in logarithmic space. The bins are ’’centered” in 2* days, 
where the bin interval is [2 l ~5, 2 l+ a] with i = 0,1, 2, 3,... 

To take into account the errors in the measurements 
we generated a random number drawn from a Gaus¬ 
sian distribution with mean and standard deviation equal 
to the measurement value and its uncertainty, respec¬ 
tively. Then, we computed the structure function using 
this mock data set and we “normalized” it removing in 
quadrature the contribution produced by the signal vari¬ 
ance estimated as the median value of the uncertainties 
in the measurements (owrf mn.')- _ 

As it is sho wn in [Schmidt et al.1 (120101! and 
iPalanaue-Delabrouille et alJ (hoilh the structure func¬ 
tion of a QSO is well described by a power-law, where r 
is measured in years, as: 

SF norm [r) = At 1 ( 6 ) 

In Figure flfll we show the power-law exponent ( 7 ) and 
amplitude (A) that resulted from the power-law fit to 
our sample. Note that A = SF{t = 1). All the ob¬ 
jects classified as variables based on the V parameter are 
marked with a white spot on Figure [T9l It is clear from 
the figure that the variability of BL AGNs is on aver¬ 
age described by larger A and 7 values compared to NL 
AGNs and GALs. The latter objects tend to be clustered 
close to 7 ~ 0.05 and A cs 0.05 in the 7 -A plane. Based 
on this we defined the region 7 >0.025 and A >0.1 in 
which the observed variability is consistent with being 
powered by the accretion disk of a super-massive black 
hole. Therefore, NL AGNs and GALs lying in this re¬ 
gion are suspected to show variability consistent with an 
AGN. 

In Figure [lH] we investigate the connection between the 
observed variability in NL AGNs and GALs and the vari¬ 
ability powered by the accretion disk of a super-massive 
black hole represented by BL AGNs. We compare ensem¬ 
ble structure functions of BL AGNs shown in blue, NL 
AGNs shown in red, and GALs shown in green when we 
apply different selection criteria to define the sample used 
to compute the ensemble structure functions. In panel 

(a) we show the ensemble structure functions when we 
consider all objects in the best observed sample; in panel 

(b) we show the ensemble structure functions when we 
consider only variable objects (V >1.3); in panel (c) we 
show the ensemble structure functions when we consider 
objects located in the BL AGN dominated region defined 
in Figure dni (7 >0.025 and A > 0 . 1 ); finally, in panel (d) 
we show the structure functions when we consider ob¬ 
jects that are located in the BL AGN dominated region 
(see Figure [TUI) and are classified as variables according 
to the V parameter. 

In panel (a) can be seen that when we consider all ob¬ 
jects in the best observed sample the ensemble structure 
functions of the NL AGNs and GALs are nearly flat, and 
very different from the AGN powered variability of BL 
AGNs that are characterized by an increase of the ampli¬ 
tude of the ensemble structure function as function of r. 
On the other hand, clear variability signatures begin to 
appear in the ensemble structure functions of NL AGNs 
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Fig. 18.— Ensemble structure functions for BL AGNs (blue), 
NL AGNs (red), and GALs (green) for different subsamples, the 
dashed-grey line corresponds to 0.05 mag the limit value from 
which any variability signal should dominate over significant noise 
on the QUEST data, (a) Structure functions for the best observed 
sample of BL AGNs (n=55), NL AGNs (n=24), and GALs (n=23). 
(b) Structure functions for the objects of the best observed sample 
classified as variables according to V >1.3, namely 44 BL AGNs, 
5 NL AGNs, 5 and GALs. (c) Structure functions for the objects 
of best ob served sample located within parameter space region of 
Figure fT9l dominated by BL AGNs (A >0.1 and 7 >0.025), namely 
48 BL AGNs, 7 NL AGNs, and 7 GALs. (d) Structure functions 
for the objects of bes t ob served sample located within parameter 
space region of Figure fl9l dominated by BL AGNs and classified as 
variables; 41 BL AGNs, 3 NL AGNs, and 2 GALs. 

and GALs when we consider variable objects (panel b) 
or objects located in the BL AGN dominated region in 
the 7 -A plane (panel c). However, only when we re¬ 
strict our selection criteria to variable objects (V >1.3) 
located in the BL AGN dominated region (7 >0.025 and 
A >0.1) we obtain ensemble structure functions for both 
NL AGNs and GALs showing variability as function r 
consistent with being powered by an AGN (i.e., similar 
to BL AGNs; see panel d). Our sample of variable ob¬ 
jects with structure function parameters consistent with 
AGN powered variability correspond to 74.5% (41) of BL 
AGNs, 12.4% (3) of NL AGNs, and 8.7% (2) of GALs. 
This implies that 3/5 and 2/5 of the variable NL AGNs 
and GAL sources selected based on the V index are now 
securely classified as variable AGNs. This is also consis¬ 
tent with the estimates of false positives derived using 
the binomial distribution. 

7.2.3. Distribution of the [3.6] — [8.0] Color 

On the right-panel of Figure[20]we show the overall dis¬ 
tribution (step-line histogram) of BL AGNs (top-panel), 
NL AGNs (middle-panel), and GALs (bottom-panel) as 
function of the [3.6] — [8.0] color. The distributions of 
variable sources with 7 >0.025 and A >0.1 (i.e., SF(t) 
similar to BL AGNs) are shown in the panels as shaded 
histograms. The grey dashed line corresponds to the 
division b etween AGN dominated and host dominated 
sources o f lBrusa et al.l d2010f) . 

As expected all BL AGNs show red colors with the ma¬ 
jority of them consistent with being dominated by the 
emission from hot dust. Variable BL AGNs do not show 
any trend with the [3.6] — [8.0] color. The [3.6] — [8.0] color 
distribution of NL AGNs shows a large spread. There is 
a group of NL AGNs with very blue colors, and a sec¬ 
ond group with redder colors peaking on the blue side of 
the BL AGN distribution, and lying close to the division 
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Fig. 19. — Amplitude (A) and power-law exponent ( 7 ) obtained 
from fitting a power law to the structure function to the best ob¬ 
served sample of BL AGNs (blue), NL AGNs (red), and GALs 
(green). We mark with a white spot the object s cla ssified as vari¬ 
ables based on the V parameter (see Section 01 . In the top- 
panel we show the distribution of the 7 for BL AGNs (blue line), 
NL AGNs (red line) and GALs (green line). In the top-panel the 
blue shaded distribution corresponds to variable BL AGNs, while 
the red and green hatched distributions correspond to variable NL 
AGNs and GALs, respectively. Similarly, on the right-panel we 
show the distributions of the amplitude (A). The black dashed 
line demarcates the region in the A —7 parameter space dominated 
by BL AGNs (namely A >0.1 and 7 >0.025), and therefore the 
region where the observed variability can be considered as being 
powered by the accretion disk of a supermassive black hole. 

be tween AGN d omin ated and host dominated sources 
of iBrusa et all (120101 ). It is important to mention that, 
although the NL AGN sample is small, among the ob¬ 
jects located on the blue side of the BL AGN color dis¬ 
tribution we find the ones with the larger excess vari¬ 
ance (see Figure El). The GALs shows a large spread 
in the [3.6] — [8.0] color distribution. However, all vari¬ 
able sources have very blue colors and are consistent with 
being host-dominated objects. 

7.2.4. Variability vs. Hardness Ratio 

The Hardness Ratio (HR) can be a measure of the X- 
ray obscuration suffered by an AGN, and it is defined as 
HR = (H-S)/(H+S), where S and H are the count rates 
in the soft (0.5 — 2.0 keV) and hard bands (2.0 — 10 keV). 
High hydrogen column densities ( Nh ) block soft X-rays, 
hence the effect of increasing column densities is that the 
X-ray spectrum becomes harder. Many studies use a Nh 
value of 10 22 cm -2 to separate be t ween obscured and 
unobscured AGNs; iMainieri et al.1 (120071 ) showed that 
90% of the sources with column densities la rger than 
10 22 cm -2 have HR >—0.3. On the other hand. lHasingerl 
(j2008l) found that using a HR value of —0.2 as a threshold 
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Fig. 20.— In the left-panel of we show the overall distribution 
(step-line histogram) of BL AGNs (top-panel), NL AGNs (middle- 
panel), and GAL ojects (bottom-panel) as function of the Hardness 
Ratio, whilst in the right-panel we show the same distributions as 
function of the [3.6] — [8.0] color. The distributions of variable 
sources are shown as shaded histograms in the panels. In the left- 
panel the grey dashed line corresponds to HR=—0.2, the division 
between unobscured (low N h) and obscured (high N h) AGNs, and 
in the right-panel the grey dashed line corresponds to the division 
between AGN dominated and host dominated sources of Brusa et 
al. (2010). 

can distinguish between obscured and unobscured AGNs. 
We will explore the relation betwee n variab ility and HR 
using the latter value, also used bv [Brusa et all (12010D . 
as a threshold. Note that the flux limit in t he soft b and 
is ~ 6 times fainter than in the hard band (IBrusa et al.1 
and that objects detected only in the soft or hard 
X-ray band have values of HR= —1 and HR= 1, respec¬ 
tively. While most of_our_GAL o bjects have HR= — 1, 
which stems from the IBrusa efTah 1 (201Cf ) criteria to clas¬ 
sify X-ray sources in the GAL class (see beginning of 
Section 7), Hasinger’s (2008) sample of galaxies shows 
typical values of —0.2< HR <0.5 (see his Figure 2). 
This can be explained as lHasingerl (|2008l ) used 2-10 keV 
detected samples only. 

In the left panel of Figure [20] we show the overall dis¬ 
tribution (step-line histogram) of BL AGNs (top-panel), 
NL AGNs (middle-panel), and GALs (bottom-panel) as a 
function of HR. The distributions of variable sources with 
7 >0.025 and A >0.1 (i.e., SF(t) similar to BL AGNs) 
are shown as shaded histograms. The grey dashed line 
corresponds to HR= —0.2, the division between unob¬ 
scured (low Nh), and obscured (high Nh) AGNs. Using 
this threshold we find that 96% ± 3% (53 ± 1) of the 
BL AGNs are unobscured, whilst 54% ± 10% (13 ± 2) 
of the NL AGNs are unobscured according to this def¬ 
inition. Of the GALs, 91% ± 6 % (21 ± 1) are charac¬ 
terized by HR< —0.2 and therefore s how a soft X-ray 
spectrum (notice that (lHasingerl 12008 ) decided to clas¬ 
sify objects without a spectroscopic AGN classification 
but with HR< —0.2 as a BL AGN, and therefore, un¬ 
der this scheme, the great majority of our GAL objects 
should be reclassified as AGN). 

Following our definition of variability consistent with 
being produced by the accretion disk of a super-massive 
black hole (V >1.3, 7 >0.025, and A >0.1) we find that 
98% ± 2% (43 ± 1) of the BL AGNs that meet these 
criteria are unobscured, 100% (3) of the NL AGNs that 
meet these criteria are obscured, and assuming that the 
X-ray emission of the GALs that meet these criteria is 
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produced by low luminosity AGN 100% (2) of them are 
unobscured. As it is expected beforehand the majority 
of the variable BL AGNs are unobscured. However, it 
is interesting that, although in the NL AGNs there is 
a large spread in the HRs and that 54% of them are 
unobscured in the best observed sample, all the variable 
NL AGNs are obscured (see FigureEOl). This result is also 
confirmed by the fact that the variable NL AGNs with 
the larger variability amplitude (oy ms >0.05 mag) are 
among obscured objects. The most likely interpretation 
is that a fraction of the accretion disk emission is leaking 
through the obscuring material. 

7.3. The XMM-COSMOS Variable Objects 

In Figure [21] we show eight light curves of the most op¬ 
tically variable sources among the XMM-COSMOS ob¬ 
jects whose observed variability is consistent with being 
produced by an AGN, including the 3 NL AGNs and 
the 2 GALs that are variables and their structure func¬ 
tion parameters are consistent with variability produced 
by an accretion disk as in the case of BL AGNs. In 
the top-right of each panel the XMM-COSMOS identi¬ 
fier number for each object is shown, while the properties 
of the objects are summarized in Table [T] All BL AGNs 
show clear variability with large amplitudes ( u rms >0.05 
mag). The NL AGNs show mid-IR colors close to the di¬ 
vision between AGN and host dominated objects and all 
have HR values consistent with gas obscuration. A likely 
interpretation for these objects is disk emission leaking 
through the obscuring material. 

Previous studies have found that between l%-4% of 
field galaxies show variabil ity that could be asso ciated 
with low-luminosity AGNs (jSaraiedini et alJfeOOOlh The 
effects of the emission produced by the accreting black 
hole on the spectrum and the images of the galaxies 
containing a low-luminosity AGNs are rarely detectable. 
For example, broad emission lines are not present in the 
galaxy spectrum even if the AGN is not obscured as the 
AGN emission is heavily diluted by the host galaxy light. 
If there is an underlying low-luminosity AGN a careful 
subtraction of two images of t he galaxy may uncove r 
the presence of a variable AGN (|Saraiedini et alJf2000D . 
Since in our survey we do not subtract a reference im¬ 
age to the rest of the images, a reliable identification of 
these low-luminosity AGNs is difficult. Given that our 
sample of galaxies is X-ray selected and some of them 
show narrow-emission lines, the probability of finding a 
low-luminosity AGN is higher than for a sample of nor¬ 
mal field galaxies. Among our sample of GALs we found 
that 22% of them show variability and 8.7% show vari¬ 
ability consistent with an accretion disk. In particular, 
XID 10569 is clearly variable (rms > 0.05 mag) and con¬ 
sistent with the variation expected for an AGN. 

8 . DISCUSSION AND CONCLUSIONS 

In this paper we present the characterization of the 
QUEST-La Silla AGN variability survey. The sur¬ 
vey is a novel effort to obtain highly sampled optical 
light curves in well studied extragalactic fields, mak¬ 
ing use of the QUEST camera at the ESO-Schmidt 
telescope at La Silla. Our observing strategy will 
provide good quality light curves for many sources 
that already have a large data set of multi-wavelength 
observations, and will increase the number of AGN 


candidates in the areas of the survey where a AGN 
classification is missing or is based only on a color- 
color selection dSchmidt et al.1 l201Ct iButler fc Bloom 


2011 

: Palanaue-Delabrouille et al. 201 it Graham et al. 

2014 

sam] 

sient 

). Additionally, the survey is producing highly 
rled light curves for all classes of interesting tran- 
s (see e.g.. Rabinowitz et al. 2012; Zinn et al.l 20141 

Cartier et al. 2015f). 

In this paper, we showed that we are obtaining good 


astrometry, with an internal precision of ~ Of.'l and with 
an overall accuracy of 0'.'2 compared to SDSS. In the 
future, the QSOs lying within our observing fields may 
be used as references to increase the internal astrometric 
precision of the survey. 

We defined the Q-band as the union of the g and r 
bands, and we showed that calibrating the Q-band as 
an equivalent to the (g + t)sdss photometric system 
yields a photometric dispersion of ~ 0.05 magnitudes. 
Furthermore, we quote, conservatively, a systematic error 
in our zero-point of 0.05 magnitudes. 

We found that the use of a linear model to fit the zero- 
point as function of the instrumental magnitude q pro¬ 
duces a better calibration of our data, increasing signifi¬ 
cantly the number of non-variable stars, and reducing the 
skewness in the Q — {g + r)sDSS residuals as function of 

Q- 

We demonstrated that we are obtaining good photome¬ 
try in the range 14.0 < Q < 19.5. Since we are observing 
a minimum of two images per tile in a given night, we 
expect to produce light curves using nightly stacked im¬ 
ages. Our preliminary results indicate that in the near 
future we will produce good quality light curves with ex¬ 
traordinary cadence of observations to a magnitude limit 
of Q ~ 20.5 (roughly r ~ 21), or deeper. 

As a way to explore the quality of the data collected 
by our survey, we studied t he optical var iability of X-ray 
XMM-COSMOS sources of iBrusa et al.l (j2010fb In this 
study, we found that the QUEST-La Silla AGN vari¬ 
ability survey is ~ 75% to 80% complete in the XMM- 
COSMOS field to magnitude limit of r ~ 20, and ~ 67% 
complete to a magnitude of r ~ 21. Therefore, the survey 
loses roughly 20% of the objects, possibly due to large 
defective areas, permanently or randomly off detectors, 
and gaps between detectors. 

Based on the variability index V (see Section m , we 
found that 80% ± 5% of the BL objects are classified as 
variable objects, while 21% ± 8% of the NL objects and 
the 22%±9% of the GAL objects are classified as variable 
objects. 

We studied the relation between a rms and the [3.6] — 
[8.0] mid-IR color, and found that objects that are redder 
in the mid-IR (possibly AGN-dominated) tend to show 
larger a rms . A possible interpretation of this result is 
that the emission of the host galaxy dilutes AGN vari¬ 
ability in objects that are more host galaxy dominated. 
We also found that overall BL objects have larger vari¬ 
ation amplitudes (e.g., higher cr rms ) than NL, and GAL 
objects. However, a significant fraction (~ 33%) of NL 
variable objects show significant variability. 

We studied the variability as function of the time 
elapsed between observations (r) in the best ob¬ 
served sample using the “normalized” structure function 
( SF norm ), i.e. after correction for measurement errors. 
We parametrized SF norm {T) as a power-law (At 7 ), and 
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TABLE 6 

Properties of Some of the Most Variable XMM-COSMOS Objects. 


Name 

Class 

z 

r mag 

p a 

1 var 

&rms 

[3.6] - [8.0] 
AB color 

Host/AGN 
dominated 15 

HR 

Obscured/ 

Unobscured c 

XID 59 

BL 

1.920 

20.50 

>0.999 

0.258 

1.41 

AGN 

-0.51 

unobscured 

XID 5607 

BL 

1.359 

20.52 

>0.999 

0.197 

1.70 

AGN 

-1.00 

unobscured 

XID 5544 

BL 

1.889 

20.21 

>0.999 

0.108 

1.54 

AGN 

-0.55 

unobscured 

XID 293 

NL 

0.445 

20.32 

>0.999 

0.137 

0.93 

AGN 

-0.02 

obscured 

XID 5047 

NL 

0.252 

20.10 

>0.999 

0.068 

0.25 

host 

1.00 

obscured 

XID 1429 

NL 

0.356 

19.77 

0.999 

0.025 

0.87 

AGN 

-0.20 

obscured 

XID 10569 

GAL 

0.348 

19.91 

>0.999 

0.067 

-1.15 

host 

-1.00 

unobscured 

XID 60406 

GAL 

0.166 

19.77 

0.978 

0.038 

-1.54 

host 

-1.00 

unobscured 


a We defined P va r = 1 — P(x 2 ) (see Section 17.21) . and is interpreted as the probability of an object to be 
variable based on its y 2 value and the number of observations. 

b We used the division between AGN or host dominated sources of IBrusa et al.l ( 12010! ') , that is based on the 
mid-IR color ([3.6] — [8.0]) of the source. 

c We follow the convention that sources with HR <—0.2 are unobscured. 
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Fig. 21.— Q-band light curves of XMM-COSMOS variable objects, 
number for each object. 

found that BL AGN tend to be located in a region with 
large amplitudes (A > 0.1) and power-law exponents 
(7 > 0.025) compared to the bulk of NL AGNs and 
GALs. We also demonstrated that the variability in vari¬ 
able NL AGNs and GALs with A > 0.1 and 7 > 0.025 
seems to be the result of an accreting super-massive black 
hole as in the case of BL AGNs and not an artifact of the 
observations. We found that 74.5% of BL AGNs, 12.4% 
of NL AGNs and 8.7% of GALs are variables and have 
A > 0.1 and 7 > 0.025. 

We confirm that the majority of the BL variable ob¬ 
jects have mid-IR colors ([3.6] — [8.0]) consistent with be¬ 
ing AGN-dominated, and hardness ratios (HRs) consis¬ 
tent with unobscured X-ray sources. On the other hand, 
we found that although 54% ±10% of the NL objects in 
the best observed sample are unobscured according to the 
HR division, 100% of the NL AGNs showing variability 
consistent with being powered by a super-massive black 
hole ( V > 1.3, 7 > 0.025 and A > 0.1) are obscured. 

This is an interesting result that requires confirmation 
from larger samples with high quality light curves. The 
time span of our QUEST-La Silla light curves will con¬ 
tinue to increase until 2016, so we expect to revisit this 
result in a future paper. 

The [3.6] — [8.0] mid-IR color distribution of NL objects 
shows a large spread, with a group of NL objects showing 
very blue colors, and a second group with redder colors 



400 600 800 1000 1200 


JD - 2,455,000.0 (days) 


In the top-right of each panel is the XMM-COSMOS identifier 


peaking close to the division bet ween AGN d omi nated 
and host dominated sources of IBrusa et al. (]2010f h but 
located on the blue side of the BL distribution (see Figure 

m- 

Although the NL AGN sample is small, we noticed 
that the objects with larger excess variance and with 
7 > 0.025 and A > 0.1 are those with relatively red mid- 
IR colors and larger values of HR (see Figure fl7l). Thus, 
a possible interpretation of variable NL objects is that 
they are obscured sources, in which the optical variabil¬ 
ity could be the result of leaked or scattered disk emis¬ 
sion. However, an increase of the sample size is definitely 
required to have conclusive interpretations. 

GALs show a large spread in their [3.6] — [8.0] color dis¬ 
tribution. However, all variable GAL sources have blue 
mid-IR colors, and thus are consistent with being host- 
dominated objects. Additionally, assuming that variable 
GALs with 7 > 0.025 and A > 0.1 contain an AGN, we 
find that these are unobscured X-ray sources according 
to the HR definition. Consequently, a plausible interpre¬ 
tation for the variable GAL objects having 7 > 0.025 
and A > 0.1 is that they are low-luminosity AGNs. 
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